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ABSTRACT 

The 6 cm formaldehyde (H2CO) maser sources in the compact H II re- 
gions NGC 7538-IRSl and G29.96-0.02 have been imaged at high resolution 
{(^heam < 50 mas). Usiug the VLBA and MERLIN, we find the angular sizes 
of the NGC 7538 masers to be ~ 10 mas (30 AU) corresponding to brightness 
temperatures ~ 10^ K. The angular sizes of the G29.96-0.02 masers are ~ 20 mas 
(130 AU) corresponding to brightness temperatures ~ 10'' K. Using the VLA, we 
detect 2 cm formaldehyde absorption from the maser regions. We detect no emis- 
sion in the 2 cm line, indicating the lack of a 2 cm maser and placing limits on the 
6 cm excitation process. We find that both NGC 7538 maser components show an 
increase in intensity on 5-10 year timescales while the G29.96-0.02 masers show 
no variability over 2 years. A search for polarization provides 3-0" upper limits 
of 1% circularly polarized and 10% linearly polarized emission in NGC 7538 and 
of 15% circularly polarized emission in G29.96-0.02. A pronounced velocity gra- 
dient of 28 kms~^ arcsecond"^ (1900 kms^^pc^^) is detected in the NGC 7538 
maser gas. 



- 2 - 



Subject headings: individual (NGC 7538) — individual (G29.96-0.02) — masers 

1. Introduction 

Astronomical observation of emission from the lio In transition at 4.83 GHz (6 cm) 
of the formaldehyde (H2CO) molecule is exceedingly rare. Only four Galactic sources have 
been observed to emit: three of these sources have been shown to be masers. The H2CO 
6 cm emission in the compact H II region NGC 7538 B was first observed by Downes 
& Wilson (1974) and was shown to be a compact, non-thermal maser by Forster et al. 
(1980). Most recently, an H2CO maser was discovered in G29.96-0.02 by Pratap, Menten 
& Snyder (1994, hereafter PMS94). H2CO masers are also observed in Sgr B2 (Whitcoak 
& Gardner 1983). The fourth H2CO emission source, Orion-KL, is known to have extended 
thermal emission {e.g. Johnston et al. 1983). Though detection of H2CO masers is expected 
to require interferometric observations which are not sensitive to the broad angular scale 
lio 111 absorption that dominates single dish measurements of the 4.83 GHz line {e.g. 
Bieging et al. 1980), dedicated interferometric surveys {e.g. Forster et al. 1985; PMS94; 
Mehringer, Goss, & Palmer 1995) have been rather unsuccessful in discovering new masers. 

There is currently no working theoretical model of the Galactic formaldehyde maser pro- 
cess, despite a history of relatively rapid understanding of astrophysical formaldehyde phe- 
nomena. Almost immediately after the discovery of interstellar formaldehyde in 6 cm absorp- 
tion by Snyder, Buhl, Zuckerman, & Palmer (1969), observations were made of "anomalous" 
6 cm absorption (Palmer, Zuckerman, Buhl, & Snyder 1969), 2 cm absorption (Evans, Che- 
ung, & Sloanakcr 1970), and emission in the 2i2 — > In (Kutner, Thaddeus, et al. 1971) and 
2ii — lio (Thaddeus et al. 1971) millimeter lines (see Fig. 1). Improved laboratory measure- 
ments (Tucker, Tomasevich, & Thaddeus 1971, 1972; Nerf 1972) and excitation calculations 
(Townes &; Cheung 1969; Thaddeus 1972; Evans 1975; Garrison et al. 1975) explained these 
phenomena shortly thereafter, but cannot explain the observed maser emission. Boland & 
de Jong (1981) proposed a continuum radiation pump mechanism to explain the NGC 7538 
maser. Though the radiative pump has been successfully apphed to the several (> 10) ex- 
tragalactic H2CO (mega)masers which have been observed {e.g. Baan et al. 1986, 1993), the 
model is not applicable to the continuum morphology discovered in NGC 7538 more recently 
{e.g. Pratap et al. 1992), nor can the model explain the Sgr B2 and G29.96-0.02 masers 
discovered subsequently. Thus, the Galactic 6 cm formaldehyde maser phenomenon remains 
unexplained nearly 30 years after its discovery. In a search for new empirical constraints on 
these mysterious objects, this paper presents an observational study of the Galactic H2CO 
masers in NGC 7538 and G29.96-0.02. 
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NGC 7538 (S158) is a well-studied H II region complex at a distance oi d ^ 3 kpc 
(Campbell & Thompson 1984). At the southern edge of the optical nebula lie three radio 
and infrared sources. The radio sources were designated A, B, & C by Martin (1973) from 2"- 
resolution 5 GHz observations. Wynn- Williams et al. (1974) detected infrared counterparts 
at 20 yUm for the radio sources; A: 1RS2, B: IRSl, and C: IRS3. The masers are located 
in front of IRSl (IRAS 23116+6111) which is partially optically thick at 6 cm. The IRSl 
H II region is thought to have a central star of type 06 V or BO. 5 II {e.g. Israel et al. 1973; 
Willner 1976) and there is a relatively large mass of cold dust associated with IRSl (but 
not IRS2 or 3) resulting in a relatively low H II region gas-to-dust mass ratio of 75 (Willner 
1976). NGC 7538-IRSl is also associated with a remarkable number of maser species besides 
H2GO (see §4.4). 

G29.96-0.02 (hereafter G29.96) is a cometary H II region at a distance of d 6.5 kpc 
(Morisset et al. 2002). The H2GO masers lie in a "hot core" about two arcseconds west of 
the cometary head. The hot core (T > 100 K; n > 10^ cm^'^) is thought to be powered by 
an embedded star, not by the cometary H II region {e.g. De Buizer et al. 2002; Pratap et al. 
1999; Cesaroni et al. 1998). The hot core also shows maser emission from H2O and CH3OH 
(Hofner & ChurchweU 1996; Walsh et al. 1998). 

All previous observational studies of H2CO masers have utilized either single dish an- 
tennas or ~ 10 km baseline interferometers {e.g. Forstcr et al. 1980, 1985; Rots et al. 1981; 
PMS94). Past observations were unable to spatially resolve the maser emission. A lower 
limit of Tb — 10^ K brightness temperature has been determined for the NGC 7538 and 
G29.96 masers (Rots et al. 1981; PMS94). The NGC 7538 H2CO source has two velocity 
components at vi,sr — —58 km s~-^ (component I) and —60 km s~^ (component II), separated 
by ~ O'.'l (Rots et al. 1981). Forster et al. (1985) noted a factor of two intensity increase in 
component I over 3 years. No polarization was detected in the NGC 7538 maser components 
to a 3-0" upper limit of 5% of the total intensity in the Stokes Q, U, and V parameters 
(Forster et al. 1985). The G29.96 masers exhibit several blended spectral features including 
velocity components at wlsr — 100 km s^^ (component I) and 102 km s^^ (component II) 
that appear coincident on the sky (PMS94). No polarization observations have been made 
of the G29.96 masers. 

Sensitive milliarcsecond-resolution observations of H2CO masers are needed. Placing 
more stringent empirical constraints on Galactic H2CO 6 cm maser emission requires the 
following observational results: (1) a determination of the angular extent of the maser emis- 
sion 6m, allowing a determination of the intrinsic sizes and brightness temperatures, (2) a 
detection of the radiatively coupled 2ii 2i2 (2 cm) line, yielding limits on inversion and 
level populations, (3) sensitive limits on polarized intensity, and (4) a characterization of the 
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observed maser intensity increase in NGC 7538. Furthermore, accurate relative positions of 
the H2CO masers with respect to other maser species could constrain the poorly-understood 
physical conditions of the maser environment. 

Compared to the Rots et al. (1981) VLA observations of the NGC 7538 H2CO masers 
which used 12 antennas recording a single polarization, the VLA is now more sensitive, using 
all 27 antennas and recording dual polarizations. Furthermore, using the sensitivity of the 
VLBA with the phased VLA, observations of the relatively weak (~ 100 mJy) H2CO masers 
at milliarcsecond angular resolution are possible. We have conducted new observations of 
the NGC 7538 and G29.96 H2CO masers at 6 cm and 2 cm using the Very Large Array 
(VLA) of the NRAO^ and at 6 cm using Very Long Baseline Array (VLBA) of the NRAO 
and the Multi-element Radio Linked Interferometry Network^ (MERLIN). 



2. Observations and Data Reduction 

We have summarized the observational parameters in Tables 1,2, and 3. Tables 1 and 2 
list the logistical parameters of the observations including bandwidth, spectral configuration, 
and calibrator information. All antennas in all observations use dual-circular polarization 
feeds. Table 3 lists the resulting sensitivities of the images. We assume throughout that 
the observed 6 cm formaldehyde line is the dominant 2 — > 2 hyperfine component which is 
separated from the 1 — > component by 1.13 km s~^ (Tucker et al. 1971; Downes et al. 
1976). 

The maser line widths, velocities, and velocity gradients were determined using the 
The Groningen Image Processing System (GIPSY) software package^. All other reduction 
and imaging was done with the Astronomical Image Processing System (AIPS) software 
package^. 

^Thc National Radio Astronomy Observatory (NRAO) is a facility of the National Science Foundation 
operated under a cooperative agreement by Associated Universities, Inc. 

^MERLIN is operated as a National Facility by the University of Manchester, Jodrell Bank Observatory, 
on behalf of Particle Physics and Astronomy Research Council (PPARC). 

■^http : //www . astro . rug . nl/~gipsy/ 

''http : //www . nrao . edu/aips/ 
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2.1. VLA 'CnB' 

The NGC 7538 and G29.96 masers were observed on 13 February 1996 at 4.83 GHz 
(6 cm) using the 27 antennas of the VLA for 20 minutes each. The NGC 7538 and G29.96 
maser positions were also observed on 13 February 1996 at 14.49 GHz (2 cm) with the VLA 
for 40 minutes each. In 'CnB' configuration the VLA basehne lengths range from 0.065 km 
to 7.7 km. 

All of the VLA 'CnB' observations were sensitive to the continuum radiation from the 
H II regions. The continuum contribution was subtracted from the maser channels using the 
AIPS task UVLSF. The 2 cm continuum radiation from G29.96 was averaged and imaged 
from 40 line-free channels over about 1 MHz of bandwidth. The RMS noise of the line-free 
2 cm continuum image of G29.96 is 3.0 mJy beam~^. 

2.2. VLA 'B' 

The NGC 7538 masers were observed on 24 July 2002 at 4.83 GHz with the VLA for 
40 minutes. In 'B' configuration the array baseline lengths range from 0.21 km to 11.4 km. 

The VLA 'B' observations were sensitive to the continuum radiation from the H 11 
regions. The continuum contribution was subtracted from the maser channels using the 
AIPS task UVLSF. The continuum radiation was averaged and imaged from 28 line-free 
channels over 84 kHz of bandwidth. The RMS noise of the line-free continuum image is 
1.6 mJy beam~^. 

2.3. MERLIN 

The NGC 7538 masers were observed on 1-3 October 2001 using the MERLIN radio 
telescope of Jodrell Bank Observatory for a total of ~ 28 hours. Six antennas were used; the 
Mark II at Jodrell Bank, the 32 m antenna at Cambridge, and the 25 m dishes at Knockin, 
Darnhall, Tabley, and Defford. The baseline lengths of MERLIN range from 11 km to 217 km. 
The phases were cahbratcd by frequent observations at 14 MHz bandwidth of 2300-1-638. 
3C84 was observed at both 14 MHz and 500 kHz bandwidth allowing bandwidth-dependent 
calibration to be transferred. 

Due to a pointing error, the MERLIN observations were centered 1.6 arcminutes from 
the maser position. Though a 2' displacement is within the primary receiving beam of the 
MERLIN antennas (FWHM ~ 10'), decoherent averaging during on-line integration ('time 
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smearing') made necessary a factor of 2.8 amplitude correction (see Thompson 1999). The 
measured positions of the masers were not affected by the error. 

2.4. VLBA+Y27 

Using the ten antennas of the VLBA and the 27-element phased VLA as an 11-station 
VLBI array, we observed the NGC 7538 and G29.96 H2CO masers on 22 September 2000 for 
~ 1.8 hours each. The amphtude scale was set by online system temperature monitoring and 

a priori gain measurements. The VLBA observations of the masers were phase-referenced to 
their respective phase calibrators. The phase referencing technique allows the maser visibility 
phases (which correspond to sky position) to be determined from those of the brighter phase- 
reference source (sec Walker 1999). The maser targets and the phasc-rcfcrcncc sources were 
alternately observed in a 2.5:1.5 minute cycle. Data from the Hancock, North Liberty, Mauna 
Kea, and St. Croix VLBA stations were not part of the final data set because the station- 
based visibility phases could not be adequately tracked and transferred due to poor signal 
to noise. The basehne lengths contributing to the final images range from 52 to 1806 km. 

3. Results 

The VLA images of the continuum radiation from the H II regions are shown in Figures 2 
and 3. The VLA 'B' array image (§2.2) of the 6 cm continuum radiation from NGC 7538 is 
shown in Figure 2 with the position of the H2CO masers indicated. The VLA 'CnB' array 
image (§2.1) of the 2 cm continuum radiation from G29.96 is shown in Figure 3 with the 
position of the H2CO masers indicated. 

3.1. NGC 7538 6 cm Emission 

3.1.1. Angular Position and Size 

We have detected H2CO masers with VLBI observations for the first time. The VLBA 
images of the NGC 7538 H2CO component I masers are shown in Figure 4. Component 
1 splits into two peaks, la and lb, in the VLBA image. The absolute position uncertainty 
of the strongest peak in the VLBA images of component 1 is 2 mas, phase-referenced to 
the calibrator J2302+6405. The positions of the emission peaks in the brightest VLBA 
channel image at i^lsr = —57.84 km s~^ are summarized in Table 4. The VLBA position for 
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component I is consistent with Rots et al. (1981) within their 100-mas-uncertainty. The la 
and lb feature peaks are separated by 14.0 ± 0.5 mas (42 AU) in the vlsr = —57.84 km s~^ 
channel image. The total extent of the velocity structure of component I in the VLBA images 
is ~ 60 mas (~ 200 AU). 

The MERLIN images show both components I and II. Figure 5 shows the channel 
images of the brightest emission from the two components. The components are separated 
by 79 ± 5 mas (240 AU), in agreement with the previous determination of 110 ± 30 mas by 
Rots et al. (1981). 

Component I shows a very compact structure from which MERLIN observations recover 
all of the emission and from which the VLBA recovers about 60% of the emission. In the 
VLBA images (Fig. 4), features la and lb have resolved angular diameters 9„i ^ 10 mas 
corresponding to brightness temperatures Tg ^ 10^ K. In contrast, component II does not 
show any compact structure on long baselines: MERLIN detects only ~ 70% of the emission 
and only the two shortest VLBA basehnes detect component II. The deconvolved size of 
component II in the MERLIN image (Fig. 5) is 64 mas. The component sizes are consistent 
with the previous upper limits of 150 mas measured by Rots et al. (1981). The angular size 
of component II corresponds to a brightness temperature of 5 x 10^ K. 

The measurement of the maser brightness temperature and of the background continuum 
radiation which the maser amplifies allows a determination of the gain of the maser. The 
measurements are related by Tb ~ T&pe""^ where Tb is the brightness temperature of the 
maser, T^g is the brightness temperature of the background continuum radiation, and r is 
the gain of the maser. Using a brightness temperature T^g — 2700 K measured from the 
VLA 'B' observations, the masers have gains r ~ —11. 

3.1.2. Velocity Gradient 

The maser fines are resolved in velocity in all of the observations except those from the 
VLA 'CnB' array. All of the fines are fit well by single gaussians. Table 5 summarizes the 
fitted and deconvolved peak intensity, velocity width at half-maximum, and line center for 
the different observations. 

The MERLIN and VLBA observations show, for the first time, a velocity gradient in 
the masing gas. In the VLBA images, the dashed arrow in Figure 4 traces the direction of 
a velocity field which extends through the la peak at a position angle of 31 ± 5° for several 
tens of milliarcseconds (several beamwidths). Figure 6 shows a position- velocity plot of the 
slice indicated in Figure 4. The plot shows a linear velocity dependence along the slice of 
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magnitude 28 kms ^ arcsecond ^ (1900 kms ^pc The MERLIN observations show a 
velocity gradient in both components that is consistent with that shown in Figure 6. 

We note that a hnear velocity gradient is consistent with the maser emission arising from 
a rigid rotating disk (or cylinder). However, the position of component II is not consistent 
with a linear extension of the component I gradient shown Figure 6. That is, component 
II does not lie on the disk which models Figure 6. Emission from a maser component at 
Vlsr = —60.19 km s"^ consistent with the 28 kms^^ arcsecond"^ velocity gradient would lie 
?a 60 mas away from component 1 at ~ 30° position angle, as marked by the '+' symbol 
in Figure 5. However, component 11 is observed to lie ^ 80 mas from component 1 at 17° 
position angle (Fig. 5). Therefore we conclude that the two maser components are not simply 
related and do not lie in a single, uniform gas cloud. 

3.1.3. Intensity 

The VLA 'B' array spectrum of the H2CO masers is shown in Figure 7. In July 2002, 
we observed component 1 to have an intensity of 1257 ± 5 mJy beam~^. The intensity on 
July 2002 is a factor of three larger than the 335 ± 20 mJy beam^^ observed in April 1983 
by Forster et al. (1985) when variabihty was first noted using the Westerbork Synthesis 
Radio Telescope with 0.36 km s~^ velocity resolution. We note an intensity increase in 
component II for the first time. In July 2002, we observe an intensity for component II of 
385 ±5 mJy beam~^ with the VLA 'B' array compared to 205 ±20 mJy beam~^ observed by 
Forster et al. (1985). A light curve of the maser component intensities is shown in Figure 8. 

3.1.4. Polarization 

At the 3-(T level, we detect no polarization in any of the observations. Our VLA data 
place a more sensitive upper limit of 1% (S-cr) on circularly polarized intensity than the 
previous hmit of 5% set by Forster et al. (1985). The MERLIN and VLBA data show no 
detectable linearly polarized intensity, corresponding to a 10% (3-0") upper limit in Stokes 
Q U, consistent with the previous limit of 5% set by Forster et al. (1985). 
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3.2. G29.96 6 cm Emission 

3.2.1. Angular Position and Size 

Table 6 summarizes the VLBA image properties of the G29.96 masers. The VLBA+Y27 
■f^LSR — 100.2 km image of G29.96 H2CO maser component I is shown in Figure 9. The 
absolute position uncertainty of the VLBA image peak is 10 mas, phase-referenced to the 
calibrator J1851+0035. The image of component II in the wlsr = 102.1 km channel 
appears coincident on the sky. The two components are separated by 1 ± 3 mas (7± 19 AU). 

The resolved sizes of the masers in the VLBA images are 9jn — 20 mas, corresponding 
to brightness temperatures Tb ~ 10^ K. Using a brightness temperature T^g = 50 K for 
the background 6 cm continuum radiation measured from the VLA 'CnB' observations after 
a beam dilution correction using the 2 cm continuum image, the G29.96 maser gains are 
r ~ -12. 



3.2.2. Line Widths 

Figure 10 shows the G29.96 H2CO maser spectra. The figure includes both the VLBA+Y27 
spectrum from the current paper and the VLA 'A' array spectrum observed by Pratap et al. 
(PMS94). Pratap et al., with a spectral resolution of 0.38 km s^^, detected "several blended 
features" whose emission is spread over ~ 5 km s^^. Remarkably, we detect two spectrally 
unresolved velocity features in the VLBA data despite a more favorable spectral resolution 
of 0.30 km s~^. Our VLA 'CnB' spectra are consistent with Pratap et al. after correction for 
absorption, which is discussed below in §3.3.2. The properties of the G29.96 hue emission 
are summarized in Table 7. 



3.2.3. Polarization 

At the 3-a level, we detect no polarization in any of the observations. The VLA and 
VLBA observations were sensitive only to Stokes / and V parameters. The VLA data place 
an upper limit of 15% (3-(j) circularly polarized intensity on component I. The VLBA data 
are also consistent with zero polarized intensity. 
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3.3. Absorption 



3.3.1. Two Centimeter 



We have detected the 2 cm transition in the direction of Galactic H2CO maser sources 
for the first time. Using the VLA 'CnB' array, we observe the 2 cm hne in absorption in 
both NGC 7538 and G29.96. All of the 2 cm hne detections have a corresponding 6 cm 
feature that is coincident in velocity. The optical depths r of all of the absorption features 
are noted in Figures 11-14. The optical depths are calculated using 



where Tl is the (negative) temperature of the line (after continuum subtraction) and ~ 
2 K is the excitation temperature of the line (Palmer et al. 1969; Thaddeus 1972). The 
derived H2CO column densities are discussed in §4. 

Figure 11 shows the 2 cm absorption at the position of the NGC 7538 mascrs. The 
2 cm minimum in the NGC 7538 observations corresponds to the position and velocity 
(f^LSR — ~60 km s~^) of (the weaker) maser component II (Fig. lib). The position of the 
2 cm minimum at the component I velocity (i'lsr — —58 km s~^, Fig. lid) lies less than two 
arcseconds away. 

Figure 12 shows the 6 cm and 2 cm spectra at the position of the G29.96 maser. No 
2 cm line is detected at the velocity of the maser but an absorption feature is seen at 
^'LSR — 107 km s~^. The 6 cm absorption evident in Figure 12a at vlsr — 107 km s~^ and 
visK — 98 km s~^ is discussed below. 

No 2 cm hne emission is detected from either NGG 7538 or G29.96 with a ?>-a brightness 
temperature sensitivity limit of 43 K. A low gain (r ~ —1) maser amplification of the 
2 cm doublet in either NGG 7538 or G29.96 would produce an emission hne of brightness 
temperature Tg > 100 K. We conclude that there is no inversion of the 2 cm line. 



Using the VLA 'GnB' array we detect 6 cm absorption in both NGG 7538 and G29.96. 
We detect 6 cm absorption in G29.96 at the position of the masers in the hot core (Fig. 12a) as 
well as against the cometary H II region. Figure 13a shows the fit to the hot core absorption 
features which arc summarized in Table 8. The maser spectrum after subtraction of the 
absorption (Fig. 13b) is consistent with the Pratap et al. (PMS94) observations. Figure 14a 




3.3.2. Six Centimeter 
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shows the spectrum at the 6 cm minimum in the VLA 'CnB' observations which corresponds 
to absorption against the cometary H II region. Figure 14b shows the 2 cm absorption 
spectrum at this position. 

We detect no 6 cm absorption at the position of the masers in NGC 7538-IRSl. Fig- 
ure 14c shows the spectrum at the 6 cm minimum of the VLA 'CnB' observations which 

corresponds to absorption against IRS2. No 2 cm absorption is detected at this position 
(Fig. 14d) which is consistent with the 2 cm observations of Mart in- Pint ado et al. (1985) 
and Evans et al. (1975) (though compare our Fig. lla,b with Fig. 9 of Martin-Pintado et 
al.). 



3.3.3. Column Densities 

The formaldehyde column densities can be calculated from the absorption spectra. Col- 
umn densities N calculated from 6 cm spectra are denoted Nq while those calculated from 
2 cm spectra are denoted A^2- The following relations were used to determine column densi- 
ties from the spectra {e.g. Wilson 1972; Wadiak et al. 1988; Osterbrock 1989) 

^ = 1.36xlO^=^J]rAT;, 

= 1.26 X 10^^22^^'" ' 

where has units cm~^, has units K, v has units km s"-*^, and where the summation 
is over the spectral channels with line radiation. The column densities are summarized in 
Table 9 grouped by cloud velocity and background continuum object. 



4. Discussion 

4.1. Absorption Clouds 

We have detected the formaldehyde distributions in NGC 7538 and G29.96 though ob- 
servations of the 2 cm and 6 cm lines in absorption against the H II regions' continuum 
radiation. The absorbing systems offer information about the ambient formaldehyde en- 
vironment. Appendix A describes the calculations which permit comparisons between the 
6 cm and 2 cm lines relating to Table 9. Section 4.2 discusses the formaldehyde calculations 
permitted by the masers lines. 
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The detection of both 6 cm and 2 cm absorption hnes toward the G29.96 hot core at 
fLSR — 107 km (Fig- 12a,b) permits a calculation of the excitation temperatures as 
described in Appendix A. The ratio of optical depths is 2.35 for which T21/T43 = 0.21. If 
we choose excitation temperatures in this ratio for the N/T^ values listed for the hot core in 
Table 9, we find comparable column densities for the 6 cm and 2 cm lines, indicating that 
the absorbing cloud is smaller than the 6 cm beam. 

The NGC 7538 mascrs' two velocity components show 2 cm absorption at both velocities 
(Fig. llb,d). Though the masers only exist in a region of diameter ~ 30 AU (Fig. 4), the 
clouds responsible for the 2 cm absorption lines have larger sizes comparable to the beam 
(~ 3000 AU). The two maser velocity components appear to come from two separate cloud 
environments. In addition to (1) the detection of the of the distinct absorbing velocities, 
evidence for two separate NGC 7538 masers includes: (2) the components have disparate 
brightness temperatures (10^ K compared to 10^ K) and maser gains (—11 compared to —6), 
(3) the components are not connected by the velocity gradient observed in the maser gas 
(Fig. 5), and (4) the masers appear to increase intensity independently (Fig. 8). We conclude 
that the masers occur in separate clouds. 



4.2. Maser Gain Calculation 

The maser amphfication process depends on the column density of the maser species. In 
this section we calculate the formaldehyde column densities from the observed maser gains. 
We use the level notation from Appendix A for our calculation. We find that a comparison 
between the column densities derived from absorption lines and maser lines indicates a 
relatively strong maser inversion. 

The gain of the maser amplification of background radiation can be expressed as {e.g. 
Elitzur 1992) 

where B21 is the Einstein S-coefficient for stimulated emission and I is the amplification 
path or gain length. Here we express r in the form of Equation 21 from Boland & de Jong 
(1981), 

T = 3.4 X 10-^3 77 x(H2CO)n(H2) ^ , 

where 77 = — (^2 — ni)/n(H2C0) such that a more negative number is a stronger inversion, 
x(H2C0) • n(H2) = n(H2C0), a;(H2C0) is the fractional abundance of ortho-formaldehyde 
with respect to the total density of hydrogen, and n(H2) is the density of hydrogen. The 
values of the parameters used in our calculations are fisted in Table 10. 
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For the values obtained from VLBA observations (r ~ —12 and Av — 0.5 km s~^, see 
Table 5), we find rj ■ a;(H2C0) • n{B.2) ■ Al = —1.8 x 10^^ cm~^. We derive limiting values 
of f] from the following arguments: (1) We consider A/ to be geometrically constrained 
to be comparable to the radius of the H II region, Rjj. For A/ = Rjj = 2.8 x 10^^ cm, 
T] ■ ^(HaCO) ■ n(H2) = -6.3 x 10"^ cm'^. (2) For NGC 7538-IRSl, van der Tak et al. 
(2000) have determined a;(H2C0) = 10~*. Adopting this value, we obtain rj ■ n{B.2) = 
—6.3 X 10^ cm~^. (3) Any non-thermal H2CO level populations must occur in the density 
range n(H2) < 6 x 10^ cm"^ (Green et al. 1978) or else the level populations 'thermalize' 
and any pumping is quenched. For instance, the Orion-KL emission lines originate in gas 
of density n(H2) ~ 10^ cm~^ which could not support non-thermal excitation (Kunter & 
Thaddues 1971; Johnston et al. 1983). For n(H2) < 6 x 10^ cm'^ we obtain r] < -0.1. The 
maser population inversion is defined not to exceed unity such that 77 > —1 which implies 
n{B.2) > 6.3 X 10^ cm^"^. The limits —l<ri< —0.1 may be compared with the typical 
inversion strength of the Boland & de Jong radiative pump model r] = —0.008. 

Alternately, we may assume that the column density = n{}l2C0) ■ Al in the maser gas 
is the same as that independently found from our absorption observations (§3.3.3). From 
Table 9, a typical value is N/T^ = 3 x lO'^^ cm'^ K'^ For = 2 K we find that 
A^ = 6 X 10~^^ cm~^. Using A^ ■ 77 = — 1.8 x 10^^ cm~^ from above we obtain rj = —0.3 which 
is consistent with the 77 value reasoned from the maser gain. 

4.3. Inversion Mechanism 

A continuum radiation excitation mechanism capable of producing the high gains ob- 
served in the 6 cm masers would also be expected to affect the 2 cm doublet level populations: 
the relatively flat Bremsstrahlung spectrum of H II region plasmas should excite the 6 cm 
and 2 cm doublets comparably (to within an order of magnitude, see also Wadiak et al. 
1988). For example, the Boland & de Jong (1981) radiative inversion model predicts a 2 cm 
maser for H II region spectra which turnover at wavelength Aq 3 cm. However, no 2 cm 
maser is present (§3.3.1) though the spectrum of NGC 7538-IRSl is observed to turnover in 
the favorable Aq ^ 3 cm range {e.g. Pratap et al. 1992; Gaume et al. 1995; Bloomer et al. 
1998; Momose et al. 2001; Akabane et al. 2001). Furthermore, since formaldehyde occurs in 
typical column densities around most H II regions (including the maser objects, e.g. van der 
Tak et al. 2000; Schoicr ct al. 2002), there is no reason on the grounds of cither continuum 
morphology or formaldehyde abundance to expect that only the few maser-b earing H II re- 
gions should support a formaldehyde inversion. Therefore, we suggest that the 6 cm maser 
inversion is due to a rare coUisional excitation, not to the readily available centimeter wave- 
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length continuum radiation which has been discussed by other authors {e.g. Pratap et al. 
1992; PMS94; Mehringer, Goss, & Palmer 1994). The remarkable paucity of formaldehyde 
masers is likely a result of selective inversion conditions such as resonant interactions with 
species exclusive to only a few H II regions. In an attempt to find unusual conditions which 
can explain the presence of coUisionally-pumped formaldehyde masers, we discuss the shock 
morphology of H II regions. 

Hill & Hollcnbach (1978) have modelled the dissociation wave and shock wave evolution 
of expanding H II regions. Our mascr environments most resemble their "Model 3" in which 
the density is n = 10'^ cm~^ and the H II region radius is R = 1.5 x 10^^ cm and the ambient 
magnetic field \B\ = 33.5 fiG of the model is closest to that of NGC 7538-lRSl (Hoare, 
private communication; Momose et al. 2001). For NGC 7538-IRSl the model correctly 
predicts the observed H2 line emissivity (~ 10~^^ erg s"^ cm"^, Bloomer et al. 1998; Hoban 
et al. 1991), shock velocity (~ 10 km s~^, e.g. Deharveng et al. 1979), and age (3 x 10^ yr, 
e.g. Martin 1973) for the H II region. 

The relevant aspect of the Hill & HoUenbach model is that the initially faster dissociation 
wave is eventually overtaken by the shock wave. The epoch at which the shock passage occurs 
is 7 X 10^ — 2 X 10^ yr, in good agreement with the age of the IRSl H II region. During 

the passage there exists a transition region about half-molecular in composition. Also, the 
region contains transient phenomena such as non-equilibrium populations of various species 
such as cool molecular gas and hot excited atoms. 

Evidence that these transition-shock coUision regions are sufficiently rare to explain the 

observed paucity of formaldehyde masers is discussed by Bcrtoldi & Drainc (1996). They 
find that for "a broad range of parameters of interest" (density, shock velocity, etc.) the 
dissociation and shock waves are merged and there is no relative evolution of the Hill & 
Hollenbach type. However, they suggest that clumps in the ambient molecular cloud are 
sites where the dissociation and shock waves are separate and where time-dependent effects 
may become important. Therefore, we suggest that formaldehyde masers occur in the small 
number of H II regions which possess (1) separately evolving dissociation and shock waves, 
(2) an age equal to the epoch of dissociation and shock wave passage, and (3) a neighboring 
molecular cloud with ambient clumps of sufficient number density and formaldehyde column 
density to cause time-dependent wave effects. 

We consider it possible that the formaldehyde maser inversion is a result of collisions 
between formaldehyde and a rare transient species within the relatively exotic conditions 
in the shock/dissociation-transition region. However, since H2CO collisions with relatively 

rare species, such as hot H atoms or excited H2 {e.g. Green et al. 1978) have not been fully 
explored, the collision cross-sections required to formulate a quantitative model are not yet 
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available. 



4.4. Other Masers in NGC 7538 

NGC 7538 is the closest of the three sources known to have H2CO masers {d = 3 kpc 
compared with G29.96-0.02 (d = 6 kpc) and Sgr B2 (d = 8 kpc)). Also, the NGC 7538-IRSl 
region contains a large number of maser species. NGC 7538 is then best suited for studies 
targeted at understanding the H2CO relationship to other maser species in close linear and 
angular proximity. For example, empirical correlations among the H2CO, OH, and H2O 
maser species have prompted searches for new H2CO masers {e.g. Forster et al. 1985). In 
the search for a working H2CO pump, an understanding of the physical conditions of IRSl 
through the observation of other masers is valuable. Studying the nearby masers that have 
relatively well-understood excitation mechanisms could considerably illuminate the H2CO 
inversion environment. 

Figure 15 shows the 22 GHz continuum image of IRSl from Gaume et al. (1995) overlaid 
with the positions of other known masers in the region. A protostellar outflow is thought to 
power the maser environment. Based on 15 GHz continuum observations, Campbell (1984) 
suggests an edge-on protostellar disk and bipolar outflow. A disk/outflow scenario is also 
supported by 22 GHz continuum and H66q; recombination line observations (Gaume et al. 
1995) which show a ~ 100 km s~^ stellar wind and apparent disk (cf. Jaffe & Martin- 
Pintado 1999). Furthermore, based on CH3OH maser observations, Minier et al. (1998) 
suggest a protostellar disk and a blueshifted southern outflow. Alternately, a shell model 
has been proposed to explain the various maser positions (Dickel et al. 1982). A cartoon 
summary of these models, including the molecular clumps that we suggest are important 
for the formaldehyde maser process, are included with the observational data in Figure 15. 
Table 11 summarizes the observational uncertainties and physical conditions of the masers. 

The ^^NHs (3,3) masers arc the only IRSl masers that could arise at the same density, 
temperature, and position as the H2CO masers (Fig. 15; Table 11). Furthermore, the velocity 
gradient observed in the ^^NHs (3,3) masers is consistent with the H2CO maser velocity 
gradient shown in Figure 5. The position of the non- metastable ^^NHg (9,6) maser in the 
IRSl region (Madden et al. 1986) is also consistent with the ^^NH3 and H2CO masers' 
positions and velocity gradient. The ^^NHs (9,6) maser species, like H2CO, is quite rare and 
poorly-understood (Elitzur 1992). Though the current position uncertainty of the NH3 (9,6) 
maser is too large to positively associate it with the H2CO, the two may be related. The 
NGC 7538 NH3 (9,6) maser has only been observed at a single epoch; however, the NH3 (9,6) 
maser in W49 is known to be variable (Wilson & Henkel 1988; Huettemeister et al. 1995). 
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Future observations of NH3 in NGC 7538 could indicate (1) an accurate absolute position 

for the NGC 7538 ^^NHg (9,6) maser near the NGC 7538 H2CO masers and/or (2) ^^NHg 
(9,6) maser intensity increase similar to the NGC 7538 H2CO masers (Fig. 8). Either (1) or 
(2) would suggest that the H2CO and NH3 masers arise in the same volume. 

The H2O masers in the region have the most disparate physical conditions from the 
H2CO (Table 11). However, recent observations (Hoare, private communication) show that 
the IRSl H2O masers appear related to the H2CO masers. MERLIN observations show 
9 H2O masers within 200 mas of the H2CO positions. Furthermore, the two H2O masers 
closest in position and velocity to H2CO maser component 11 (fLSR — 60.1 km s~^) have 
tripled in intensity in 10 years. It is interesting to note that the H2CO masers have also 
tripled in intensity over 10 years. Since it is unlikely that the H2O and H2CO masers exist 
in the same volume, perhaps the two species are related by a shock interaction with a clump 
at vlsr — 60.1 km s~^. Indeed, the shock-driven pumping scheme for H2O masers (Elitzur, 
Hollcnbach, & McKee 1989) includes molecular clumps and may be compatible with the 
H2CO pumping mechanism we suggest using the Hill & HoUenbach (1978) and Bertoldi & 
Draine (1996) models. 



5. Conclusions 

We have detected and imaged the unusual H2CO masers with VLBI baselines for the 
first time showing milliarcsecond structure (tens of AUs in linear scale) and brightness tem- 
peratures in excess of 10* K. MERLIN images of the NGC 7538 H2CO masers show the 
separation between the two velocity components to be 79 ± 5 mas and the components ap- 
pear to arise in separate clouds. We also measure a 1900 kms^^pc^^ velocity gradient in 
the NGC 7538 maser gas. Both NGC 7538 maser components have been shown to increase 
intensity on 5-10 year timescales despite remaining constant in intensity for the first 10 years 
after discovery. We do not detect variability in the G29.96 masers on a two year timescale. 

We detect the 2 cm formaldehyde line in absorption in both NGC 7538 and G29.96. 
We detect no 2 cm line emission from any of the maser regions from which we conlcude 
that there is no inversion of the 2 cm levels. From the good agreement between the column 
densities derived from the maser lines and the absorption lines, we find that the inversion of 
the masers is 77 ~ —0.3. We suggest that the strong but rare formaldehyde inversion process 
involves a coUisional pumping interaction. We posit that the environment which generates 
these pump conditions results from unusual shock dynamics around the H II regions. 

I.M.H. is supported by the NRAO pre-doctoral researcher program. We thank D. S. 
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Shepherd for helpful comments on an earlier draft. We thank M. Hoare for communicating 
the MERLIN observations of the NGC 7538 H2O masers. 

A. Optical Depth Calculations 

The optical depth r for thermal absorption is proportional to 

\Avy \gL guj 

where /i is the permanent dipole moment of the molecule, ly is the transition frequency, Av 
is the thermal Doppler line width, S is the line strength, and n and g are, respectively, 
the number of absorbers and the statistical weights of the upper and lower levels of the 
transition. For the 6 cm and 2 cm transitions of formaldehyde Sqctu ~ 3/2, S2cm ~ 

5/6, 

(7 = 3 for both levels in the 6 cm doublet, and = 5 for both levels in the 2 cm doublet (see 
Townes & Schawlow, 1955). 

We may analyze our data using a four level system for which we use the following 
notation: In = 1, lio = 2, 2i2 = 3, and 2ii = 4. We shall present our results as constraints 
on three excitation temperatures; T21 between the 6 cm doublet levels, T43 between the 2 cm 
doublet levels, and T31 between the J = 2 and J = 1 doublets. We examine only those 
absorbing systems for which we have detected both 6 cm and 2 cm absorption. Also, we 
assume that the absorbing clouds are smaller than the telescope beam while the extent of 
the continuum radiation is not. 

When both the 6 cm and 2 cm line have been observed we may compare the optical 
depths of the lines 

< r2e,„ > ^ Qecrn 5 2 3 N, {l - e-^'-^^Vfe^^^) 

< Teem > ^2cra 6 3 5 iVi (1 - e-^^^^i/^^^O 

t^iS N3 T21 



^2cm 3 Z/21 Ni T43 

where it is assumed that hu -C kT, where fl is the solid angle of the telescope beam over 
which the absorption measurement is averaged, and where 

Ni gi 

can be substituted since the total populations between the 6 cm and 2 cm doublets are 
approximately thermally distributed. Substituting 1^43/1^21 — 3, we have the relation 

< '^2cm > ^ ^6cm 5 721 ^-husi/kTsi 

< Teem. > ^2cm 3 T43 
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For calculations we will use ^6cm/^2cm — 9.35, T31 = 20 K, and 1/31 — 140.84 GHz. 
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Fig. 1. — Energy level diagram of the lowest rotational states in the ground vibrational 
state of formaldehyde (H2CO). The level notation is Jk'K+- The K~ = 1 level splitting is 
exaggerated for presentation. The solid arrows show the radiative electric dipole {AK' = 0, 
AK~^ = ±1) transitions between the states. Note that i^~-even and i^~-odd ladders are 
not radiatively connected; the K~ — ladder is shown for perspective but does not enter 
our calculations. The K~ — 2 and K~ — 3 ladders which begin at ~ 60 K and 
E/k 120 K, respectively, are not important for this paper. The numbers in itahcs label 
the transition energies in the K~ — 1 ladder in frequency units of gigaHertz. The insets 
magnify the 6- and 2-centimeter transitions of interest. 
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Fig. 2.— The 4.8 GHz continuum radiation observed from NGC 7538 with the VLA 'B' 
array. The contours are -4, 4, 8, 16, 24, and 32 times the noise RMS 1.6 mJy bcam~^. The 
beam, which is plotted in the lower left hand corner, is 1.3x1.0 arcseconds at a position 
angle of —10°. The crossed lines mark the position of the H2CO masers. 
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Fig. 3. The 14.5 GHz continuum radiation observed from G29.96-0.02 with the VLA 'CnB' 
array. The contours arc -5, 5, 10, 20, 30, 40, 50, and 60 times the noise RMS 3.0 mJy beam~^ 
(no contours at -5 appear). The beam is 2.0x1.4 arcseconds at a position angle of 37°. The 
crossed hues mark the position of the H2CO masers. 
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Fig. 4. — Six adjacent channel images from the VLBA+Y27 observations of NGC 7538 
maser component I. The contours are -5, 5, 7, 9, 11.5, 14, and 16 times the noise RMS 
12.8 mJy beam~^ (no contours at -5 appear). The beam, which is plotted in the lower right 
hand corner, is 13x6 mas at a position angle of 37°. The dashed arrow is at position angle 
31°. marking the direction of the linear velocity field fitted to Figure 6. The 'X' marks the 
position corresponding to the bright la peak in the wlsr = —57.84 km s~^ channel image. 
Properties of the emission features are summarized in Table 4. 
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Fig. 5. — MERLIN images of the component I and component II H2CO maser line emission 
from NGC 7538. The contour levels arc 25, 50, 80, 100, and 130 times the image RMS noise 
9 mJy beam~^. The beam for the MERLIN images is 53 x 50 mas at position angle 39°. 
The dashed arrow represents the direction of the velocity gradient apparent in the VLBA 
data (Figure 4). The '+' symbol marks the position at which a fLSR = —60.19 km s~^ maser 
would lie if consistent with the extrapolated VLBA velocity gradient (Figs. 4 & 6). 
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Fig. 6 — Position-velocity plot from the VLBA+Y27 observations of NGC 7538 feature la 
showing the plane marked by the dashed arrow in Figure 4. The contour levels are the same 
as in Figure 4. The emission is fit by a linear velocity field of 1900 kms~^ pc~^ as shown by 
the dotted arrow. 
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Fig. 7.— VLA 'B' array Stokes I spectrum of the NGC 7538 H2CO masers from July 2002 
after continuum subtraction. 
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Fig. 8. — Intensity liistory of tlie NGC 7538 H2CO masers. From tlie current paper, only 
the July 2002 VLA 'B' and February 1996 VLA 'CnB' observations are plotted. 
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Fig. 9. — VLBA+Y27 image of the vlsr = 100.2 km channel (component I) from the 
G29.96-0.02 observations. The contours are -3, 2, 3, 4, 5, 6, 7, 8, and 8.6 times the noise 
RMS 7 mJy beam~^. The image of component II (i'lsr — 102 km s~^) appears at the same 
sky position. The beam is 23 x 17 mas at a position angle of —17°. The images properties 
of the masers are summarized in Table 6. 
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Fig. 10. — Spectra of the H2CO maser in G29.96-0.02 as observed with the VLA 'A' array 
by Pratap et al. (PMS94) (soUd hne) and the VLBA+Y27 (dashed hne). For the Pratap et 
ah VLA observations the spectral resolution is 0.38 km s^^ and the beam size is 550 mas 
(0.02 pc). For the VLBA observations the spectral resolution is 0.30 km s~^ and the beam 
size is 20 mas (0.0006 pc). 
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Fig. 11. — VLA 'CnB' spectra at 6 cm and 2 cm of the NGC 7538 H2CO masers after 
continuum subtraction. Spectra (a) and (b) are both taken at the position of the maser 
maximum which is also the position of deepest 2 cm absorption (aj2ooo — 23 13 45.37, 
<^j2ooo — 61 28 10.4). Spectra (c) and (d) are both taken at the position of the vlsr — 
—58 km s~^ (component I) 2 cm minimum (q;j2ooo — 23 13 45.71, 5j2ooo — 61 28 10.7). The 
beam size for the 6 cm observations is 5" (0.07 pc), for the 2 cm observations is 2" (0.03 pc). 
The T values are the peak optical depths of the respective lines. 
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Fig. 12.— VLA 'CnB' spectra at 6 cm and 2 cm of the G29.96-0.02 H2CO mascrs after 
continuum subtraction. Spectra (a) and (b) are taken at the same sky position (q;j2ooo = 
18 46 03.73, (^j2ooo — ~02 39 22.5): the location of the 6 cm maser maximum. The absorption 
in (a) Bit vlsr — 98 km s'^ and i'lsr — 107 km s"^ is subtracted in Figure 13b. The beam 
size for the 6 cm observations is 5" (0.16 pc), for the 2 cm observations is 2" (0.06 pc). The 
T values are the peak optical depths of the respective lines. 



-34- 



_ ' 1 ' 
- (a) 

X = 0.62 


1 ' 
\; x= 1.86 


1 ' 


1 1 1 

- J 

^^^^^^^ 
, ^ 1 


-l 1 

r 

1 


1 1 — 

1 



140 



120 



100 



80 



Fig. 13. — VLA 'CnB' 6 cm spectra for G29.96-0.02. (a) shows the observed spectrum of 
Figure 12a (sohd histogram) and the fitted absorption profiles (dotted curve), (b) shows the 
spectrum in (a) after subtraction of the fitted absorption summarized in Table 8. The r 
values are the peak optical depths of the respective lines. 
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Fig. 14. — Comparison of VLA 'CnB' 6 cm and 2 cm absorption spectra. The position of 
(a) Sz (b) is the location of the 6 cm minimum in the G29.96-0.02 observations (6 cm beam 
is 5" (0.16 pc), 2 cm beam is 2" (0.06 pc)). The position of (c) (d) is the location of the 
6 cm minimum in the NGC 7538 observations (6 cm beam is 5" (0.07 pc), 2 cm beam is 2" 
(0.03 pc)). The r values are the peak optical depths of the respective lines. 
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Fig. 15. — {lef't) The 60-mas-resolution 22 GHz continuum VLA image of NGC 7538 B 
(IRSl) from Gaume et al. (1995) overlaid with the positions of known masers. The contours 
are -4.75, 4.75, 9.5, 14.25, 19, 23.75, 28.5, 33.25, 38, and 42.25 times the image RMS noise 
0.4 mJy beam-^ The ^^NHg, OH (4765 MHz), and OH (1720 MHz) masers are not plotted 
individually, but are shown as regions. The absolute position uncertainty of the continuum 
image is 100 mas. Due to large positional uncertainties, the ^^NHs and OH (6035 MHz) 
masers which may lie in this region are not shown. Positional uncertainties, physical con- 
ditions, and references for the masers are in Table 11. {right) Cartoon summary of the 
proposed models for IRSl. 



Table 1. Observed Rest Frequencies 



Transition Frequency 



lio ^ 111 (6 cm) 4.829657 GHz 
2ii ^ 2i2 (2 cm) 14.488475 GHz 



Table 2. Observational Parameters 



Array 


Band 




bandwidth 


bandwidth nchan 


amp/fringe/band cal'r'' 


phase calibrator 






(km s-l) 


(km s-l) 


(kHz) 




Source 


S (Jy) 


Source 


s{Jy) 


a$ (mas)'' 




NGC 7538 


VLA 'CnB' 


2cm 


-58.0 


32.0 


1562.5 


63 


3C286 


3.40 


B2229+695 


1.1 


<2 


VLA 'CnB' 


6cm 


-58.0 


95.0 


1562.5 


127 


3C286 


7.47 


B2229+695 


1.1 


<2 


VLA 'B' 


6cm 


-59.0 


7.0 


98.0 


63 


3C48 


5.48 


B2229-I-695 


1.1 


<2 


MERLIN 


6cm 


-59.0 


30.0 


500.0 


256 


3C286 


7.47 


623004-638 


1.0 


100 


VLBA+Y27 


6cm 


-59.0 


30.0 


500.0 


128 


B1749-f096 


2.2 


J2302-f6405 


0.25 


<2 


G29.96-0.02 


VLA 'CnB' 


2cm 


100.0 


32.0 


1562.5 


63 


3C286 


3.40 


B1749+096 


2.1 


<2 


VLA 'CnB' 


6cm 


100.0 


95.0 


1562.5 


127 


3C286 


7.47 


B1821-I-107 


1.3 


<2 


VLBA+Y27 


6cm 


100.3 


30.0 


500.0 


128 


B1749-I-096 


2.2 


J1851-f0035 


0.4 


10 



''The VLA and MERLIN use this continuum source to set the amplitude scale and to characterize the bandpass. The VLBA uses this 
source to search for the clock-like station delays which produce fringes and to characterize the bandpass. 

''The absolute position uncertainty of the phase calibration source as listed in the VLA Calibrator Manual at 
http : //www . aoc . nrao . edu/~ gtaylor/calib . html. 
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Table 3. Image Parameters 



Array 


Band 


Date 


spec. res. A/rms 


beam 


beam P. A. 


n a 
t7LAS 








(km s~^) (mJy beam~^) 


(as) 


(deg) 


(as) 




NGC 7538-IRSl 


VLA 'CnB' 


2cm 


13 Feb 1996 


0.60 6 


1.9x1.5 


53 


90 


VLA 'CnB' 


6 cm 


13 Feb 1996 


0.80 6 


5.7x4.5 


38 


120 


VLA 'B' 


6 cm 


24 Jul 2002 


0.20 5 


1.3x1.0 


-11 


36 


MERLIN 


6 cm 


02 Oct 2001 


0.14 9 


0.053x0.050 


38 


1.2 


VLBA+Y27 


6cm 


23 Sep 2000 


0.30 13 


0.013x0.006 


37 


0.25 


G29.96-0.02 


VLA 'CnB' 


2cm 


13 Feb 1996 


0.60 6 


2.1x1.5 


22 


90 


VLA 'CnB' 


6cm 


13 Feb 1996 


0.80 6 


6.4x4.6 


11 


120 


VLBA+Y27 


6cm 


23 Sep 2000 


0.30 7 


0.023x0.017 


-17 


0.25 



^largest angular scale to which the observations are sensitive 
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Table 4. NGC 7538 ^;lsr = -57.84 km VLBA Image Features 



Feature 


R.A.^ 


Dec.^ 


I 


^max 


Tb 




(h m s) 


(° ' ") 


(mJy beam~^) 


(mas) 


(10« K) 


la 


23 13 45.3550(3) 


61 28 10.442(2) 


216(13) 


9(3) 


1.6(2) 


lb 


45.3568(3) 


10.447(2) 


153(13) 


11(4) 


1.4(2) 



Coordinates in J2000 epoch 
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Table 5. NGC 7538 Maser Line Profiles 



Feature 


(mJy beam~^) (km s~^) 


Av 
(km s~^) 




VLA 'B' 


I 

II 


1268(20) -58.05(1) 
372(20) -60.22(8) 


0.99(6) 
0.88(8) 


MERLIN 


I 

II 


1310(35) -57.90(4) 
270(35) -60.19(6) 


0.82(5) 
0.64(8) 


VLBA+Y27 


la 
lb 

II 


216(13) -57.90(2) 
153(13) -57.82(2) 

a 


0.37(9) 
0.41(9) 



^feature not imaged 
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Table 6. G29.96-0.02 VLBA Image Features 



Feature 


R.A.^ 


Dec.^ 


I 


^max 


Tb 




(h m s) 


(° ' ") 


(mJy beam"^) 


(mas) 


(10^ K) 


I 


18 46 03.7911(7) 


-02 39 22.221(10) 


61(8) 


18(7) 


9.2(12) 


II 


03.7910(7) 


22.221(10) 


23(8) 


20(9) 


3.5(12) 



Coordinates in J2000 epoch 
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Table 7. G29.96-0.02 Maser Line Profiles 



Feature 


I vlsr 
(mJy beam^^) (km s~^) 


Av 
(km s""^) 




VLA 'CnB' 


I (& II) 


64(26) 100.43(5) 


2.42(64) 


VLBA+Y27 


I 

II 


70(18) 100.24(4) 
19(7) 102.0(2) 


0.288(87) 

a 



^feature not spectrally resolved 
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Table 8. G29.96-0.02 Hot Core 6 cm Absorption 



^^LSR 


Av 


/ 


(km s~^) 


(km s~^) 


(mJy beam~^) 


98.20(120) 


3.4(15) 


-58(22) 


107.11(90) 


a 


a 



^feature not spectrally resolved 
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Table 9. Absorption Column Densities 




(km s"-*-) 


N/T^ (from 6cm) N/T^ (from 2cm) 
(10^3 cni-2 K-i) (10^^ cm-2 R-^) 


Figure 




NGC 7538-IRSl 


-57.6 
-58.1 
-60.1 


0.27 
2.42 
3.36 


lib 
Ud 
lib 


NGC 7538-IRS2 


-57.3 


3.60 


14c 


G29.96-0.02 Hot Core 


107.1 
98.2 


2.17 7.79 
18.99 


13a, 12b 
13a 


G29.96-0.02 Cometary H II Region 


100.8 
98.5 


0.71 

0.61 0.78 


14a 
14a,b 



no absorption detected 
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Table 10. Maser Gain Parameters 



Parameter 


Symbol 


Value 


line width 


Av 


0.5 km 


gain length 


Al 


2.8 X 10^^ cm 


H2CO fraction 


^(HaCO) 


10-8 


gas density 


n(H2) 


6 X 10^-5 cm-3 
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Table 11. NGC 7538-IRSl Maser Parameters 



Maser 


Model Conditions 


Fig. 15 Position 


Species 


T 

-'■ kin 


logio n 




Uncertainty 




(K) 


(cm~^) 


Ref. 


(mas) 


Ref. 


H2CO 


20 


< 5 


1 


2 


§3.1.1 


i^NHs (3,3) 


> 50 


4-5 


2,3 


100 


10 


OH (1720 MHz) 


90 


> 5 


4 


100 


11 


OH (4765 MHz) 


~ 90 


> 4 


5,12 


500 


12 


OH (6035 MHz) 


> 170 


> 7 


6 


1800 


13 


CH3OH 


> 50 


7-8 


7 


30 


14 


H2O 


400 


7-9 


8 


15 


15,16 


i^NHa (9,6) 


~ 300 


4-10 


9,17 


3000 


17 



Note. — References: (1) Boland & de Jong 1981, (2) 
Walmsley & Ungerechts 1983, (3) Flower, Offer, & Schilke 
1990, (4) Guibert, Elitzur, & Nguyen-Q-Rieu 1978, (5) 
Elitzur 1976, (6) Pavlakis & Kylafis 2000, (7) Sutton et al. 
2001, (8) Elitzur, HoUenbach, & McKee 1989, (9) Mauers- 
berger, Henkel, & Wilson 1987, (10) Gaume et al. 1991, (11) 
Forstcr ct al. 1982, (12) Palmer, Gardner, k Whiteoak 1984, 
(13) Baudry ct al. 1997, (14) Minicr, Booth, & Conway 2000, 
(15) Kameya et al. 1990, (16) Hoare, priv. comm., (17) Mad- 
den et al. 1986 



